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Abstract 

We develop a systematic algorithm for constructing an jV-fold supersymmetric system from a 
given vector space invariant under one of the supercharges. Applying this algorithm to spaces of 
monomials, we construct a new multi-parameter family of AA-fold supersymmetric models, which 
shall be referred to as "type C" . We investigate various aspects of these type C models in detail. It 
I/-) ' turns out that in certain cases these systems exhibit a novel phenomenon, namely, partial breaking 

of TV-fold supersymmetry. 
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I. INTRODUCTION 



A significant progress in scientific research has often been achieved through the unifi- 
cation of seemingly unrelated concepts. Recently, three different theoretical developments 
were unified in the framework of M-fold supersymmetry 0, 0] , namely: i) isospectral trans- 
formations, traced back to the work of Darboux in the late nineteenth century j^] (see also 
Ref. |H and references therein), and their higher-derivative generalizations first formulated 
in Ref. [H[ , ii) quasi-exact solvability in one-dimensional quantum mechanical systems [f| 
(see also Ref. [7( and references therein), and iii) a particular class of nonlinear superalge- 
bras. The characteristic feature of A/"-fold supersymmetry, which distinguishes it from other 
nonlinear extensions of ordinary supersymmetry such as parasupersymmetry 0, 0, ^| and 
fractional supersymmetry [Hj . is the fact that anticommutators of fermionic operators are 
polynomials of degree (at most) M in bosonic operators. Usually, A/"-fold supercharges are 
represented by A/th-order linear differential operators (see references in Ref. 

The unification mentioned in the previous paragraph comes about schematically as fol- 
lows. An A/"-fold supersymmetric quantum system always yields, by definition, a pair of 
isospectral Hamiltonians H^. These operators are automatically quasi-solvable, since they 
preserve the kernel of (the bosonic part of) the respective supercharges Q ± . Finally, the 
anticommutator {Q + ,Q~} is a polynomial of degree (at most) M in the superHamiltonian 
H, which accounts for the polynomial character of the superalgebra closed by H and Q^. 
Although these ideas are conceptually simple, the direct construction of AAfold supersym- 
metric models becomes rather unwieldy for large A/", since intertwining relations with respect 
to a higher-order differential operator are quite complicated. Indeed, type A A/"-fold super- 
symmetry 0, Q] is virtually the only known class for which most of the aspects mentioned 
above are well understood for arbitrary Af '. 

One of the characteristic features of type A A/"-fold supersymmetric models is that, after 
a suitable gauge transformation, their supercharges leave invariant the space of polynomials 
in one variable of degree less than A/". In our previous paper [15[ we obtained a new family 
of A/"-fold supersymmetric models, the so-called type B, by considering what is probably the 
simplest deformation of the type A A/"-fold supercharge. Although the construction presented 
in this reference was rather ad hoc and of a purely analytic nature, the resulting type B 
supercharges (after an appropriate gauge transformation) preserve a finite-dimensional linear 
space of monomial type. These results strongly suggest that it can be of great advantage 
to base the construction of A/"-fold supersymmetric models not on a specific form of the 
supercharge, but rather on the finite-dimensional linear space invariant under one of the 
supercharges, say Q~ . 

In this article we show that this idea is indeed feasible, by developing a systematic 
algorithm for constructing an A/"-fold supersymmetric system starting from the knowledge 
of the M- dimensional linear space of functions left invariant by a suitable gauge transform 
of the supercharge Q~ . One of the main advantages of this method is that it completely 
bypasses the "hard" calculation of the intertwining relations, which makes it ideally suited 
when M is not fixed a priori. We then use our algorithm to derive the A/"-fold supersymmetric 
systems arising when one chooses as the starting point in the construction what is perhaps 
the simplest possibility, namely a space of monomials. According to Post and Turbiner 
there exist essentially three inequivalent finite-dimensional monomial spaces preserved by 
a nonzero second-order linear differential operator. It turns out that two of them lead to 
the already known A/"-fold supersymmetric models of types A and B. The remaining one 
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yields a new type of supersymmetry, referred to as type C in what follows, which is one 
of the main contributions of this paper. The type C A/"-fold supercharge is in a sense the 
most general supercharge of arbitrary order A/", since it formally reduces to all previously 
known instances, namely, the type A and B supercharges. We fully classify the type C 
models, finding in particular explicit formulas for their potentials, and further analyze the 
main properties of these models, such as shape-invariance, associated polynomial families of 
Bender-Dunne type, and the resulting AAfold superalgebras. 

The article is organized as follows. In the next section, after reviewing briefly the con- 
cepts of AAfold supersymmetry and quasi-solvability, as well as their mutual relationship, 
we outline the general procedure for constructing an A/"-fold supersymmetric system starting 
from a given finite-dimensional linear space of functions. The connection with the construc- 
tion in Refs. [3 EH of type A and B models is also briefly discussed. In Section IIHI we 
apply the general method to the simple case in which the given linear space of functions is of 
monomial type, obtaining a new multi-parameter family of AAfold supersymmetry, namely 
type C. We investigate the general properties of type C models, such as the shape invariance 
between the partner Hamiltonians, the structure of the solvable sectors, and the symmetry 
transformations which preserve the potential form. Using the invariance under these sym- 
metry transformations, we completely classify the type C models in Section HVl We find that 
there are essentially four inequivalent nontrivial types of potentials. The normalizability of 
the solvable sectors is also briefly examined in this section. In Section El we study the poly- 
nomial families of Bender-Dunne type associated with the AAfold supersymmetric models 
of type C. We prove that to each type C model one can associate two polynomial families 
which, in contrast to their type A counterparts, are always weakly orthogonal. They also 
exhibit a novel feature, namely their dependence on two integer parameters. It is also shown 
that the polynomial part of the type C AAfold superalgebra can be expressed as a product 
of two critical polynomials belonging to each of the above families. A few examples are 
also exhibited explicitly at the end of this section. Finally, in Section IVT1 we discuss several 
general aspects of type C models, such as the partial breaking of AAfold supersymmetry, 
and briefly discuss future developments suggested by the present work. 

II. AA-FOLD SUPERSYMMETRY AND QUASI-SOLVABILITY 

In this section we shall review the concept of AAfold supersymmetry in one-dimensional 
quantum mechanics, making special emphasis on its connections with the recently intro- 
duced notion of quasi-solvability. We shall then outline a general algorithmic procedure for 
constructing AAfold supersymmetric models starting from an Af- dimensional linear space 
invariant under the action of a second-order linear differential operator. 

Let q denote a bosonic coordinate, and let ip an d ^ be fermionic coordinates satisfying 



{^^} = {f I f} = ) 
Given a monic A/th-order linear differential operator 



{V,#} = 1. 



(II.l) 



Af-l 




(II.2) 



we introduce the AA-fold supercharges by 



(II.3) 
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where the operators Pjy are defined by 



P M = P N , Pb = {-Vfl*„. (II.4) 

The superscript * in the latter equation denotes the transposed operator defined by A 1 = 
(A*)*, where the star denotes complex conjugation. If, as is usually the case, all the coeffi- 
cients Wk in Eq. pi. 2)1 are real, then Pjj obviously coincides with Pjj. The nilpotency of the 
fermionic variables ip an d ip' implies that 

{QJr,QJr} = {Qtr>Qjsr} = °- ( IL5 ) 

We define a superHamiltonian H by 

H = H'tp^ + H+ipty , (II.6) 

where 

H ± = -\dl + V ± (q) (H.7) 

is a pair of scalar Hamiltonians. An M-fold super symmetric model is a triple (H, Qj^-, Qjj) 
such that the supercharges commute with H, namely 

[QJ,H]=0. (II.8) 



Equation (|II.8Jl is equivalent to the following intertwining relations between the component 

AT 



supercharges Pjb and Hamiltonians H ± 



P M H- - H + P u = 0, P^H+-H-P^ = 0. (11.9) 

Since the Hamiltonians are both symmetric under transposition, each one of the relations 
(III.9J) actually follows from the other one by transposition. 

Following Ref. 3|> we shall say that a differential operator T is weakly quasi-solvable 



with respect to an TVth-order differential operator P/y of the form (jll.2|) if it leaves ker P/v 
invariant. If, in addition, kerP^ can be explicitly computed we shall simply say that T is 
quasi-solvable. We shall also use the term quasi- exactly solvable to refer to a quasi-solvable 
operator whose finite-dimensional invariant space is a subspace of the Hilbert space on which 
the operator is naturally defined. 1 A particular class of quasi-solvable operators is that of 
Lie-algebraic operators, which are polynomials in the generators of a finite-dimensional Lie 
algebra of first-order differential operators preserving a known finite-dimensional linear space 

00,1110. 



± 



From Eqs. (jll.9|) it immediately follows that each one of the component Hamiltonians H 
leaves invariant the kernel Vjy- of the corresponding operator P^. Hence both H~ and H + are 
(weakly) quasi-solvable with respect to the operators P^ and P^, respectively. Remarkably, 
the converse is also true. Indeed from a scalar Hamiltonian H weakly quasi-solvable with 



1 Unfortunately, the definition of the term "quasi-exactly solvable" in the literature is far from uniform. In 
this paper we have adopted the terminology recently proposed by one of the authors in Refs. 0,0. 



4 



respect to an A/th-order linear differential operator Pjy of the form (|II.2|) one can construct 
an jV-fold supersymmetric model by taking Pj^ as in (|II.4|) and 



H = H 



H + = H + w M -M 



(11.10) 



where the dot denotes derivative with respect to q. Thus there is a one-to-one correspondence 
between A/"-fold supersymmetric models and weakly quasi-solvable scalar Hamiltonians. 
From the above remarks it follows that, given an TV-dimensional linear space 



V^= (<Pi(q),...,w(Q)) 



(11.11) 



and a scalar Hamiltonian H leaving Vj^ invariant, one can construct an A/"-fold supersym- 
metric model using Eqs. (jll.3jl - (jll.4j) and (jll.lOj) and taking as Pv the unique A/th-order 
linear differential operator (|II.2|) annihilating V/j-, namely 



m) 1 



1 d q 



q 

(AO 



(11.12) 



where WUpj , ■■ ■ ,(f \r) is the Wronskian of the functions (px,...,<pjy spanning W 
(cf. Refs. |20lI21||). The latter scheme for constructing A/"-fold supersymmetric models suffers 
from two major drawbacks, namely: i) an arbitrary A/"-dimensional linear space (jll.llj) is 
not preserved, in general, by any scalar Hamiltonian, and ii) even when this is the case, it 
may be difficult to find explicitly a Hamiltonian leaving (jll.llj) invariant. 

To overcome these drawbacks, let us slightly generalize the above construction by con- 
sidering an M- dimensional linear space 

Vu = {<pi(z),...,<PM(z)) (11.13) 

and a scalar second-order linear differential operator (not necessarily a Hamiltonian) 

-H = A{z)d 2 z +B{z)d z + C{z) (11.14) 

leaving W invariant. From these two ingredients one can construct an A/"-fold supersym- 
metric model as follows. 
Let 



Pm = g(z) 



V k=0 



w k {z)d] 



(11.15) 



denote the most general A/th-order linear differential operator with kernel VV, where the 
function g(z) is for the time being undetermined. We shall first construct another second- 
order linear differential operator of the form 



H + =H-SC 



(11.16) 



satisfying the intertwining relation 



P N H -H + P, 



(11.17) 
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To this end, note that the left-hand side of (|II.17|) is in general a linear differential operator 
of order Af + 1. Equating to zero the coefficients of d^f +l and d^f in this operator we obtain 
the following two equations for the functions g and 5C: 



i 



5C 



J\f(J\f - 2) A" - 



A 
2A 



+ J\f ( B' - 



BA' 
~2A 



- A'wm-i - 2Aw N _ 



i • 



(11.18) 



(11.19) 



where the prime denotes derivative with respect to z. When Eqs. (|II.18|) - (|II.19j) are satisfied, 
the l.h.s. of (|II.17|) is a linear differential operator of order at most M — 1 annihilating the 
A/"-dimensional space VV, and hence it vanishes identically. 

The last step in our construction consists in applying a change of variable 



z = z{q) 

and a gauge transformation with gauge factor e~ w ^ z \ under which 



H 



± 



z=z(q) 



(11.20) 



f 11.21^ 



to simultaneously take H~~ = H and H + to Schrodinger form (|II.7jl . Note that this is 
certainly possible, since (by construction) H and H + differ by a scalar function only. The 
appropriate change of variable and gauge transformation are determined by |l9i . |22| 



z z = 2A{z) 



W 



( 



.4' 



2A V 2 



B 



(II.22a) 
(II.22b) 



The potentials V ± (q) are related to the coefficients of the differential operators if ± as follows: 



V ± (q) = 
where C~ = C, C H 



- C± + Ta 
AA 



B' 



AA" + - A' 2 + 2{AB' — A'B) 



(11.23) 



C~ + SC. More explicitly, 



V ± (q) = -C + 
1 



A'wm- 



Aw 



z=z(q) 



A" 



+ 



AA 



Q 2 + (A/" 2 



A" 



± 



A 

MQ' (HQ w X -i \ 

■i- — + \Oa + -z- ]A 



(11.24) 



z=z(q) 



where, following Ref. |15l |. we have set 



1 



Q = B + -{N-2)A'. 



(11.25) 
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From the above construction it immediately follows that the system (|H.3|) - (jH.4|) and (|II.6|) 
with H^f given by (HL2TJ) and 



(11.26) 



is jV-fold supersymmetric. Indeed, the first intertwining relation (|II.9|) follows by applying 
the gauge transformation and change of variable (jll.21|) to the relation (jII.17|) . Note also 
that (jll.26|) and the definition of Pj^ imply that the kernel of the operator Pjj- = Pjsf is given 
by (jnTEJ), with 

<Pi(q) = e" W(z) fr(z)\ z=z(q) , \<i<N. (11.27) 

Likewise, the invariance of the space VV under the Hamiltonian H~ is an immediate conse- 
quence of the invariance of VV under H and Eqs. (jll.21|) and (jll.27|) . It is important to note 
that Eq. (ITL26J) for P N is compatible with Eq. (ITL21 . Indeed, from Eqs. (iHTsl) and (lll.22al) 
it follows that g is proportional to . Taking 

9{*) = %)"U(*) (H.28) 

and using Eqs. (jIL15|) and ()II.26|) . it immediately follows that P^ is indeed of the form 
flTOl. It i s also straightforward to check, using Eqs. (ITHl) . (IILT5|), (ITLT9|) . (IIL23|) . (ITL"26j) . 
and pi. 28)1 . that the partner Hamiltonians H ± are related by (|II.10jl . 

At this point we shall briefly summarize the results obtained so far. Given an M- 
dimensional linear space (III. 13 J) and a second-order linear differential operator H =H~ 
(III.14J) leaving it invariant, one can construct an A/"-fold supersymmetric model (H, Q^f) 



through the following algorithmic steps: 

i. Compute the change of variables z(q) using Eq. (jll.22a|) . 

ii. Construct the operator 



N — Z 



M-l 

d?+J2Mz)d k z ) (IL29) 



k=0 



whose kernel is the linear space VV- 

iii. Compute the operator H + pi,16j) using Eq. (|II.19j) . 

iv. The component Hamiltonians H ± are obtained from the gauged Hamiltonians H ± by 
applying the gauge transformation (jll.21|) . with W(z) given by ()H.22b|) . 

v. Likewise, the operator P^ determining the A/"-fold supercharges through 
Eqs. flll.3jl - pi.4j) is obtained from P^ via the same gauge transformation, 
cf. Eq. (JTISHJ). 

In practice, the component Hamiltonians H ± (|II.7|) are computed using Eq. ()IL24j) for the 
potentials V ± (q). Note also that the second step in the previous construction is algorithmic, 
since the operator in parenthesis in Eq. (jll.29)) can be computed using a formula analogous 
to (jII.12|) . In most cases, however, the operator Pjy is easily found by inspection. 

To make the above construction completely symmetric with respect to the partner Hamil- 
tonians H~ and H + , and to facilitate comparison of the above results with previous work 
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[1J, [l5| , we next introduce two additional differential operators P5- and H + as follows. In 
the first place, note that from the second intertwining relation in Eq. (jll.9|) it follows that 
H + leaves invariant the kernel of the supercharge 

P£ = (-lW = (-l) M e™P^e-™. (11.30) 

Using the identity 2 

(dzY=(-d q ) =-d q \ = -zd z \ 
\z J z z 

and Eq. (jTL29jl we can write 

{-lY r P* N = z 1 -"Pfri»- 1 , (11.31) 

where 

Af-l 



P^ = z^ldf +J2(-^- k d k z w k ) . (11.32) 
\ k=o J 

From Eqs. (pT30jl and fTT^Tj) we have 

^ = e" W+ ^e w+ , (11.33) 

where the function W + is given by 

W + = -W + (A/" - 1) In \z\ . (11.34) 
The invariance of ker P£ under H + and Eq. (jll.33|) imply that the operator 

H+ = e w+ H + e' w+ (11.35) 

leaves the linear space 

V£ = kerJft (11.36) 
invariant. Setting W~ = W we can express the partner Hamiltonians as 

H± = e- w± H±e w± , (11.37) 

where the gauged Hamiltonians if ± leave invariant the kernel of the "gauged" supercharges 
Ptf respectively given by Eqs. ()H.29|) and ()H.32j) . The "physical" supercharges Pjjj are 
related to the gauged supercharges by the equations 

P% = e~ w± P±e w± . (11.38) 
In order to express W ± in a symmetric way, we introduce the functions 

W(q)= l -{w-(q)-W + (q)) (11.39) 



Note that the transposition has been defined in terms of the variable q, cf. Eq. (|II.4 
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and 

E(q) = Z M . (11.40) 

z{q) 

From Eq. (jll.22a|) . its immediate consequence 

z — A' , 

and Eq. (JII.22bj) it is straightforward to derive the relation 

W = -Q. (11.41) 

z 

We then have 



= J E(q)dqT J W(q) dq = - 1) ln\2A(z) \ ± f£fifite, (11.42) 

cf . . The connection between the gauged Hamiltonians H + and H + follows easily from 
Eqs. (jnTq) and ([Tl35jl . namely 

= e~ 2 / w (9) d 9#+ e 2 / w(9)d<? . (11.43) 

Using Eqs. (lH~T9l) . (I11.22al) . (jTL25j), and ([11.41 jl . it is immediate to obtain the following 
explicit expression forH + : 

-H + = A8 2 z -(q + I(jV- 2)A'J d z + C + (M - l)Q' - A'wjf-i - 2Aw' JV _ 1 . (11.44) 

Combining this equation with Eq. (jII.14j) we obtain the following unified formula for the 
gauged Hamiltonians H ± : 

-H ± = Ad 2 z - (±Q + ^{N-2)A^j d z + C 

+ 1(1 ± 1) {{M-l)Q' - AW-i - 2Aw' Ar _ 1 ) . (11.45) 

This general formula includes as particular cases the gau ged Hamiltonians of the type A and 
type B models introduced respectively in Refs. [l^ and |l5j. Indeed, in the type A models 
we have Wj^f-i = and 

C = i(jV - - 2) A" l)Q' + R , 

where R is a constant. Using these relations in Eq. (jll.45[) we can easily reproduce Eqs. (3.41) 
and (3.50b), (3.52), (3.55) of Ref. [lij for the gauged Hamiltonians of type A (note that 
the coordinate h and the function P in the latter reference correspond to our z and A, 
respectively). Likewise, type B models satisfy w^f-i = —1/z and 

C=±-{N-l){N-2)A"-^-\{N-l)Q'--^ r + R, 

12 2z 2 J\l z 

which immediately lead to Eq. (3.35) of Ref. [l5j for the gauged Hamiltonians of type B. 



9 



III. A NEW FAMILY OF AA-FOLD SUPERSYMMETRIC SYSTEMS 



In this section we shall apply the previous results to the construction of a new multi- 
parameter family of TV-fold supersymmetric systems for arbitrary Af. The key idea in this 
respect is to choose appropriately the A/"-dimensional linear space ()H.13j) . in such a way that 
the linear space of second-order linear differential operators leaving it invariant is nontrivial 
and can be explicitly computed. 

To this end, we shall consider monomial spaces of the form 

V M = (z x \... 7 z x «), (III.l) 

where the exponents Aj are real numbers. All spaces (jllllj) left invariant by a nonzero 
second-order linear differential operator have been classified by Post and Turbiner 161 . up 
to changes of variables and gauge transformations of the form 

ip(z) i-> $(z) = z a ip(z p ) , a,/3eR. (111.2) 

For Af > 4 the above classification consists of three equivalence classes, represented by the 
following "canonical forms" : 

A) (M,...,^- 1 ) (III.3) 

B) <l,z,...,^- 2 ,^) (III.4) 

C) (l,z,...,z Ari - 1 ,z x ,z x+ \...,z x+Af2 - 1 ), Af = Af 1 +Af 2 . (III.5) 

In the third canonical form Af\ and Af 2 are positive integers, and A is a real number different 
from -Af 2 , -Af 2 + 1, . . . ,Aft (A / -Af 2 -l,Af 1 + l when Kf x = 1 or Af 2 = 1 to prevent the 
canonical form C from reducing to B). Due to the freedom in performing changes of variables 
and gauge transformations (jIII.2|) we can also assume, without loss of generality, that 

A > , Afx>N 2 . (III.6) 

The A/"-fold supersymmetric models constructed from the canonical form A following the 
procedure described in the previous section are nothing but the type A systems introduced 
in Ref. [13]. Similarly, the models obtained from the second canonical form are the type B 
systems recently constructed by the authors . In this section we shall therefore derive and 
completely classify all the A/"-fold supersymmetric models associated to the last canonical 
form (|III.5J) . which from now on we shall term type C models for short. 

We shall now proceed to the construction of the type C models by following the algorith- 
mic steps outlined in Section |H] In order to implement this algorithm, we must know at 
least one nonzero second-order linear differential operator leaving the space (|III.5|) invari- 
ant. Due to the extremely simple nature of the latter space, it is actually straightforward 
to determine the whole linear space of linear differential operators of order not greater than 
two preserving it. Indeed, using the techniques described in Ref. j3] it is readily found 
that the latter space is spanned by the constant multiplication operator 1 and the following 
operators: 

(HI.7a) 
(III.7b) 

Af 2 + 1) (HI. 7c) 

(III.7d) 



Jo- 


= d z (zd z 


-A) 


Joo 


= z*dl 




J+o 


= z(zd z - 


■ Mi + 1) (zd z - A 


Jo 


= zd z . 
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The most general linear second-order differential operator leaving the space (jIII.5|) invariant 
is thus given by 

-H = a% J - + a 2 Joo + a 3 J +0 + Mo + c , (III.8) 

where the coefficients a%, bo, and Co are real constants. More explicitly, the coefficients A(z), 
B(z) and C(z) of —H (cf. (111.14)) ) are given by 

A(z) = a 3 z 3 + a 2 z 2 + a\z (III. 9a) 

B(z) = -(TV + A - 3) a 3 z 2 + b z + (1 - A) a x (III.9b) 

C(z) = (Mi - 1)(/V 2 + A - 1) a 3 z + c . (III.9c) 

By Eq. (ITL251) we have 

Q = -(JV- 2A) a 3 z 2 + b lZ + ~(Af - 2A) a x , (111.10) 

with 

b 1 = b + (Af-2)a 2 . (Hill) 

From Eqs. (|II.22a|) and (|IH.9aj) . the change of variable z(g) is determined in this case by the 
differential equation 

2(a 3 z 3 + a 2 z 2 + a x z). (111.12) 



z 



2 



The type C space (IIII.5)) decomposes in a natural way as the direct sum of two spaces of 
type A, namely 

(1, z, . . . , z^- 1 ) ®z x (l,z,..., z^- 1 ) . (111.13) 

It is important to observe that both subspaces in the latter sum are separately invariant under 
all of the operators (IIII.7|) . and hence under the gauged Hamiltonian H. In particular, a 
result of Turbiner |23J] implies that H and z x H z~ x are Lie-algebraic operators with respect 
to the standard realization of the algebra sl(2) with generators 

J_ = d z , J = zd z , J + = z 2 d z -nz (111.14) 

and cohomology parameter n = M\ — 1 and n = M 2 — 1, respectively (this property can 
also be checked directly using (|III.7j) ). However, an arbitrary polynomial in the operators 
flIII.14|) will not preserve, in general, both type A spaces in flIII.13|) . For this reason, the 
number of independent first- and second-order operators preserving Vj^ is reduced from 8 
for type A to 4 for type C. In particular, the set of type C gauged Hamiltonians does not 
include all gauged Hamiltonians of type A. In Section El we will discuss this phenomenon 
from a different viewpoint, namely the breakdown of an underlying symmetry. It should also 
be noted in this respect that the gauged Hamiltonian of a type B model is not, in general, 
a polynomial in the sl(2) generators (|HI.14|) . 

The operators J _, Joo, an d Jo obviously leave invariant the spaces (1, z, . . . , z^ 1 ^ 1 ) 
and z x (l,z, . . . ,r /v ' 2_1 ) for arbitrary positive-integer values of M\ and M 2 . It follows from 
Eq. f)III.8|) that the gauged Hamiltonian H will also leave invariant the latter spaces for all 
Af±,Af 2 G N provided that the coefficient a 3 vanishes. When this is the case, Eq. (111.37)) 
implies that the Hamiltonian H~ preserves the two infinite ascending sequences of spaces 

V K = e- w -^z(qf-^ x (l, z(q), . . . , z^" 1 ) ; i = 1, 2 , A^GN. 
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Hence for 03 = the type C component Hamiltonian H is solvable in Turbiner's sense 
z4, zfc> see also [13 E3) We will see shortly that in this case the other Hamiltonian H + 
is simultaneously solvable, as for type A models. 

The A/"th-order linear differential operator P^ of the form (|II.29j) having as kernel the 
type C space (|III.5J) is easily found to be 

Pm = i N [d z + ^^V- (ni.i5) 

In particular, the function wj^-i is given in this case by 

Wtf-i = . (III. 16) 

z 

To compute the operator P^ determining the supercharges Q^, let us introduce the 
function 

F(q) = 44 , (HI.17) 
related to E by the identity (see Eq. (jll.4f)|t ) 

F = EF-F 2 . (111.18) 

Making repeated use of the equality 

(d z + -) = i" 1 ^ + aF)i^ fc = z-* _1 (a, + fc£) (III. 19) 
we immediately obtain 

AA-i A/i-i 

*V = J] (3, + (.A/i - A)F - iS) J] (9, - iS) , (111.20) 

where the products of operators are ordered according to the following definition: 

h 

J~J = ^h^h-i ■ ■ ■ ^4*0 • 

i=io 

Using the identity 

e~ w 9 9 e w = d q + W = d q + ^(M- l)E(q) + W(q) (111.21) 

and Eq. 26)1 we finally obtain the following explicit formula for the operator P/y of type 
C: 

AT-1 



= J] (d q + W + (Af 1 -\)F + ^(Af-l-2i)E^J 

A/i-i / 1 \ 

J| (d q + W+-(J\f—l- 2i)E j . (111.22) 



X 

i=0 
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It is clear from the previous expression that the supercharge ()III.22|) reduces to its type A 
counterpart for A = M\ = Af — 1, and to the type B one for A = A/i + 1 = M . We can 
thus formally regard the type C supercharge as a deformation of those of types A and B 
depending on the two parameters A/i and A. 

The pair of type C potentials V ± can be expressed in terms of the functions E, F and 
W using the following identities, which are easily derived from Eqs. (|II.22aj) . (|II.40J) . (|II.41J) . 
rfTTTTBl). and (jTTTTTjl : 



= A4(M - X)EF , 2Aw' Af _ 1 = -N 2 {Nx - A)F 2 , (111.23) 

(111.24) 



2A = z 2 F\ ^ = E 2 , A" = E + E 2 , 



O ON 

2/i = W2 > ^A = ~ EW > Q' = -(W + EW). 

Substituting the above formulas into Eq. (|II.24|) and using the relation piI.18jl . we finally 
obtain 3 



2 q 2 4 24 V J 



2 2 



where R is a constant. 

From Eqs. ()H.32j) and (jIII.15j) it follows that the operator P£ is given in this case by 

p^=^(d M -^^y a . (IH.25) 

Its kernel Vj^ is easily computed (cf. Eq. (|II.36jl ). with the result 

= ^ 2 (1, z, . . . , z U \ . . . , ^ +A/ ' 2 - 1 ), (111.26) 

where 

A = A/i - AT 2 - A . (111.27) 

The space can obviously be transformed into the type C canonical form (jIII.5j) by the 
gauge transformation ip(z) i— > ^(z) = 2~- /v ' 2 -?/'(z). Since A is negative if A > J\f\ — A2, strictly 
speaking it is necessary in this case to perform an additional gauge transformation and 
change of variable ()HI.2j) to take the space ()IH.26j) into the canonical form (|III.5j) - (|III.6J) . 
It is preferable, however, not to enforce the first restriction A > in Eq. (jHI.fij) in the 
sequel (cf. the discussion following Eq. (jIII.41jl ). the positivity of A being immaterial in 
other respects. Hence from now on we shall only impose on A the restriction needed to 
ensure that the type C space (jIII.5|) does not reduce to the type A or B canonical forms, 
namely 

AgK\{-A/ 2 ,-A/' 2 + 1,...,M}, (111.28) 

with A ^ -Af 2 - l,A/i + 1 if A/i = 1 or Af 2 = 1. It is easily seen from Eq. 1)111.2711 that A is 
restricted in exactly the same way as A. 



Due to the relation (jIII.18fl . this expression is not unique. 
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From Eq. (|III.26|) it follows that the gauged Hamiltonians H + [ai, 6 1( c , A] and 

a>i, bi, Cq, A] are related by the gauge transformation 



H + [a h 61, c , A] = z-^F [tii, 61, c , A] z 



-A/2 



(111.29) 



for a suitable choice of the parameters aj, 61, and Co- The latter equality implies, by the 
uniqueness (up to an additive constant) of the physical Hamiltonian associated to a given 
gauged Hamiltonian, the important relation 



H + [ai,bi,co, A] = H [a^ 61, c , A] . 



(111.30) 



In other words, the type C A/"-fold supersymmetric models we shall obtain are guaranteed 
to be (formally) shape invariant j2^|. Explicitly, the parameters (cij,&i,Co) are given by, 



1,2,3, 



61 = -61 + 2A/" 2 a 2 , 

c {i = c + (N 1 -l)(b 1 -N 2 a 2 ). 



(III.31a) 

(m.31b) 
(III.31c) 



Since a 3 = implies 03 = and vice versa, it follows that H~ and H + are always simulta- 
neously solvable. 

We shall next examine the explicit forms of the subspaces Vjj- preserved by the type C 
Hamiltonians H ^ . We first note from Eqs. (jIII.9a|) and (jIII.10|) that the second term of the 
last expression for W ± in Eq. (|II.42J) is given by 



Q(z) , J\f -2\ fh U-2X 

2Mz) dz = —r~ ln N + U " 4 



A(z) 



dz 



Introducing the new parameters 

1 



a 



-A 



a+=7V 2 -A/'i + A + ^ = i- A. 



(111.32) 



(111.33) 



we obtain the following expression for the gauge factors of the type C models: 



a. 



1 ± 1 



W 2 ln\z\ + 

2 2 211 



Af-1 



In 



A(z) 



M-2X 

± a 2 

2 4 



dz, (111.34) 



where an irrelevant constant term has been dropped. By Eq. (|II.37|) . the subspaces Vjy 
invariant under the Hamiltonians H ± , which provide the algebraically computable wave 
functions (not taking their normalizability into account), are obviously given by 



From Eqs. (JIlLBj) . (I111.26j) and ()111.34j) we finally have 



(111.35) 



V± = e- U± z^ a± 



^- 1 )©e" w± ^ 1 - ± )(l,,,...,/- 1 ) 



2=2(9) 



(111.36) 
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where the new gauge factors e u± are defined by 

TV - 2A 



A(z) 




exp 


z 





b_i 
2 



-a 2 



A(z 



dz 



(111.37) 



The last step in our construction is the computation of the component Hamiltonians H ± 
using Eqs. (jll.7|) and (jll.24j) . This is, technically speaking, the most delicate step, since it 
involves the explicit evaluation of the elliptic integral 



[2(a 3 z 3 + a 2 z 2 + a\Z 



-1/2 



dz = ±(q - q ) 



(111.38) 



needed to compute the change of variable z = z(q), cf. Eq. piI.12j) . 

The value of the integral (|III.38J) . and hence the corresponding change of variable, depends 
on the position of the roots of the polynomial A(z) in the complex plane. We should therefore 
classify the polynomial A(z) into (real) canonical forms according to the position of its roots, 
using changes of variables and gauge transformations 

i;(z)^^(z) = ii(z)^(z)\ z=as) (HI.39) 

that preserve the form of the type C space f)III.5j) . This task is hindered by the fact that, 
contrary to what happens in the analogous classification of the one-dimensional Lie-algebraic 
and type A A/"-fold supersymmetric Hamiltonians 0, 0, 2^1 , the space (jIII.5j) is not in- 
variant under translations z = z + Zq. Fortunately, however, this space is invariant under 
dilations 

i/>{z) i — ► ${z) = ip{az) , a e R , (111.40) 
and is form-invariant under special projective transformations 

t/j(z) \-+$(z) = z s ^{z" 1 ) , s = Mi -I, JV 2 + A-1. (111.41) 

Since, however, neither projective transformation pil.41|) preserves both conditions (|III.6|) . 
we must drop one of these conditions if we insist on using projective transformations to bring 
A(z) into canonical form. From now on we shall assume that only the second condition 
Mi > N 2 in Eq. (jlll.fij) holds, which entails the choice s = Mi - 1 in Eq. (|HI.41jl . 

The transformations f)III.40j) — f)III.41() map the gauged Hamiltonian H (jII.14j) into the 
operator H respectively given by 



-H = -H 



Q- 



A(az) 



dz< 



+ - B(az 



a 



d_ 

di 



G{az) 



and 
-H 



-z s Ez 



A(z 



^± + s2 ) 2 + B(z- 1 )( -z 2d 



dz^'V 1 )(-^ + sz)+C(z-i). 

In particular, the polynomial A(z) transforms under dilations and special projective trans- 
formations respectively as 



and 



A{z) 1 — ^ A{z) 



A{z) 



1 



a- 



A(az) 



(111.42) 



A{z) = z A A{z~ l ) . (111.43) 

With the help of the transformations pil.42jl - pil.43jl it is readily shown that A(z) can be 
cast into one of the canonical forms listed in Table UJ The discussion following Eq. pil,12j) 
implies that the type C models corresponding to the first two canonical forms in Table HI or 
to the third one with m = 0, are not only quasi-solvable but also solvable. 
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TABLE I: Canonical forms for the polynomial A(z) (|III.9a|) . In this Table v > 0, < m < 1 and 
m' = 1 — m (m ^ 1 in Case 3 and m ^ 0, 1 in Case 4 to avoid duplications). 



1) 2z 

2) ±\vz 2 

3) ±2vz{l - z)(l-mz) 

4) 2zvz(l — z)(m' + mz) 

5) ii/z(z 2 + 2(l-2m)z + l) 



IV. CLASSIFICATION OF THE TYPE C MODELS 

We shall now explicitly compute the type C models associated to each of the canonical 
forms in Table HJ Note that, by Eq. (|III.12|) . a rescaling of the coefficients di, b±, cq by an 
overall nonzero constant factor v has the following effect on the change of variable z(q): 

z(q; Vdi, z/fei, vc ) = z(y/uq; a h b u c ) . (IV.l) 

From this equation and Eqs. (fTT24"j) . (lH~25j) . (|lX^ - (|lL42)l . and (I111.17j) we easily obtain 
the identities 

E(q; va h vb u uc ) = y/uE(y/uq; a i} h, c ) (IV.2a) 

F(q; vcii, vb\, uc ) = \puF{^Jvq\ d^ b x , c ) (IV.2b) 

W(q; vcii, ubx, ucq) = y/vW(y/vq, di, b h c ) (IV.2c) 

W ± (q; va h ub u uc ) = W ± (^q; d t , b u c ) (IV.2d) 

V ± (q; vdi, vh, vcq) = v V ± (^/uq; a i} h, c ) . (IV.2e) 

We shall therefore set u — 1 in the canonical forms 2)-5), the models corresponding to an 
arbitrary value of v following easily from Eqs. (jlV.lJ) and (jIV.2j) . It should also be obvious 
from Eq. (jIII.38|) that the change of variable z(q), and hence the functions E, F, W and the 
potentials V ± determining each model, are defined up to the transformation q \— > ±(q — q ), 
where go £ ^ is a constant. We shall make use of this observation to simplify the expressions 
for E, F, W, and V ± . 

Case 1. A(z) = 2z. 

Change of variable: z = q 2 . 

Supercharge: 



E=K F=-, W= -h q+ l^^l. (,V,) 
q q 2 2q 



Potentials: 



1,2„2 , ^(O^-I) N 

2q 2 + 4 

Here, and in what follows, Vq denotes an arbitrary constant 



V ± = ±6? q 2 + - ^ ; =F + V . (IV.4) 
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Solvable sectors: 

V^=e^ 2 q a± (l,q 2 ,...,q 2 ^) , (IV .5a) 
V^=^q l -^{l A \...^^) , (IV.Sb) 

where the linear spaces Vjy. are defined in Eq. (jIII.36|) . This case corresponds to a solvable 
model. The potentials (|IV.4J) are singular at the origin, and thus their Hamiltonians may be 
naturally defined on, e.g., the half-line S = (0, oo). In this case, we see from Eqs. (IIV.5J) that 
the normalizability of the solvable sectors depends on the value of or* 1 (or A) and the sign of 
b\. The finiteness of the L 2 norm in the solvable sectors yields the following conditions: 4 

V^cL 2 (S) 61 > , A > Ni - H 2 - 1 , (IV.6a) 

V^CI 2 (S) 61 > , A < M - N 2 + 1 , (IV.6b) 

V M C L 2 (S) 61 < , A < 1 , (IV.6c) 

V N2 CL 2 (S) 61 < , A > -1 . (IV.6d) 

Case 2a. A(z) = z 2 /2. 
Change of variable: z = e 9 . 
Supercharge: 

E = F = 1 , W — —b\ . (IV. 7) 

Potentials: 

V ± = V . (IV. 8) 

Thus both potentials are equal and trivial in this case. 
Case 2b. A(z) = -z 2 /2. 

The formulas for the supercharge and the potentials for this case can be easily deduced from 
those of the preceding one by applying Eqs. (jIV.2|) with v — — 1. 

Case 3a. A(z) = 2z(l - z)(l - mz). 

Change of variable: z = sn 2 q. 

Here, and in the following cases, the Jacobian elliptic functions have modulus k = y/m. 
Strictly speaking, if < m < 1 the above change of variable is only valid in one of the two 
regions in which A(z) is positive, namely the interval < z < 1. In the second region of 
positivity 1/m < z the change of variable is z — l/(m sn 2 g). However, since the projective 
transformation w = l/(mz) leaves A invariant and maps the interval (0, 1) into the half-line 
(1/m, 00), we need not consider the second change of variable. 
Supercharge: 

3msn 4 g — 2(1 + m) sn 2 q + 1 cng dng 

E = , F = 2 , 

sn g cn g dn g sn g 

(IV.9) 

(1 — Af — 2a )(1 + m sn g) — bi sn g 
2 sng cng dng 



4 For the solvable sectors to be included in a Hilbert space on which the Hamiltonians H ± are self-adjoint, 
we need of course to impose a stronger restriction on the parameters coming from the boundary condition 
at the endpoint q — 0, which is usually taken as lim g ^o+ <7 ^VK?) = 0- ^ ee also the discussion in Section 

eh 
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Potentials: 

„± "I,,, ,, , a ± (a ± -l) m'P ± (p ± -l) 

■ ^, : ') t f f + t (rv.io) 



Parameters: 



Solvable sectors: 



2 dn 2 g 2 



0* = ^(1 - AO ± (6i + (1 + m){Af+ 2a - 1)) . (IV.ll) 



= (sngf (cngf (dngf^l, sn 2 g, . . . , (sn g) 2(M - 1} ) , (IV.12a) 
V^ 2 = (sng) 1 - Q± (cng) /3± (dng) /3T <l,sn 2 g,...,(sng) 2(A/ ' 2 - 1) ) . (IV.12b) 



It is worth mentioning that in this case the potentials V + and V are related by a complex 
translation (up to a constant term), namely 

V + (q) = V-(q + IK') + ^L- (/T - /?+) . (IV.13) 

In the latter equation K' = K(w!) is the complete elliptic integral of the first kind, defined 
by 

r /2 alt 

K (m) = / . (IV.14) 

Jo v 1 — m sin 2 1 

In general, the potentials (jiV.lOj) have real singularities at the points q = nK (n £ Z), and 
thus their Hamiltonians may be naturally defined on, e.g., 5 = (0, if). In this case, we see 
from Eqs. (jIV.12j) that the normalizability of the solvable sectors depends on the value of 
and /3 ± . The restrictions on the parameters coming from the finiteness of the L 2 norm 
in the solvable sectors are now given by 

VjcI 2 (S) ^ a±> -\> (3±> -\' (IV - 15a) 
VjcI 2 (S) <=► ot<\, /? ± >-^- (IV.15b) 

As previously discussed, when m = the polynomial A(z) is of second degree, and hence 
the Hamiltonians are solvable. The formulas for the supercharge, the potentials and 
the solvable sectors are obtained from Eqs. (jiV.9|) - (jIV.12j) by setting m — 0, w! — 1, and 
(sn q, cn q, dn q) = (sin q, cos g, 1 ) . 

Case 3b. A(z) = -2z{l - z)(l - mz). 

As in Case 2b, the formulas for this case can follow from those of the preceding one using 
Eqs. flIV.2|) with v = — 1. The following well-known identities j2?J may be of help in this 

case: 

. sn(g;m') 1 dn(g; m') 

sn(ig; mj = l — , cn(ig; m) = — , dn(ig; m) - 



cn(g; m') ' 1 cn(g; m') ' cn(g; m') 
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The value m = yields again solvable models, whose supercharge, potentials and solvable 
sectors follow from Eqs. (jIV.9|) - (|IV.12j) by setting m — 0, ml — 1, and (snig, cnig, dnig) = 
(i sinh q, cosh q, 1). This case deserves further discussion, since the resulting potentials 



2 sinh^ q 2 cosh 2 q 2 ' v 2 



are now of hyperbolic type, and hence are not periodic on M. The solvable sectors are given 
by 

= (sinh q) a± (cosh q) p± (l,sinh 2 g, . . . , (sinhg) 2 ^ 1 ^ 1 )) , (IV. 17a) 

V% 2 = (sinh g) 1 -^ (cosh q) p± (l,sinh 2 g, . . . , (sinh g) 2 ^ 2 ^) . (IV. 17b) 

The potentials fjlV.161) are singular only at the origin, and thus their Hamiltonians may be 
naturally defined on, e.g., S = (0, oo). The finiteness of the L 2 norm in the solvable sectors 
leads to the following restrictions: 

C L 2 (S) ^ -- < a* < -p* - 2M + 2 , (IV.18a) 

VjcL 2 (S) (3 ± + 2Af 2 -l<a ± < ^ (IV.18b) 

Note that the above inequalities cannot be satisfied unless (3^ < —2Mi + 5/2, % = 1,2. 
Case 4. A(z) — 2z{A — z)(m' + mz). 
Change of variable: z = cn 2 g. 

Again, the above change of variable is only valid in the interval < z < 1. In the second 
region of positivity of A(z), namely the half-line z < —m'/m, the correct change of variable 
is z = —m'/(m cn 2 g). As before, we shall restrict ourselves to the interval < z < 1, 
since the projective transformation w = —ml / (mz) maps this interval into the half-line 
(— oo, —m'/m) and leaves A invariant. 

Supercharge: 



3msn 4 g — 2(1 + m) sn 2 g + 1 sng dng 

E = : , F = — 2 



W 



sng cng dng cng 
' 1 — M — 2a~)(m cn 4 q — m') +b\ cn 2 g 



(IV.19) 



2 sng cng dng 
Potentials: 

. m _ _ . 2 ^±(^±-1) m' 0^(0^-1) 

-T^T^ + T (M-M)a- + ^ + Vo. (IV.20) 
2 dn g 2 z 

Parameters: 

= i(l - AT) ± ~ (6i + (1 - 2m)(M + 2a- - 1)) . (IV.21) 
As in the previous case, the scalar potentials are related by a complex translation: 

V+(q) = V~(q + K + iK') . (IV.22) 
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Solvable sectors: 

V± = (cng) Q± (sng) /3± (dng) /3T <l,cn 2 g,...,(cng) 2 ( A/ ' 1 - 1 )), (IV.23a) 
y± = (cng) 1 - a± (sng) /3± (dng) /3T <l,cn 2 g,...,(cng) 2 ( A, ' 2 - 1 )). (IV.23b) 

The potentials are again singular at integer multiples of K, so that their Hamiltonians H ± 
can be taken as defined on, e.g., S = (0,K). The conditions ensuring the finiteness of the 
L 2 norm in the solvable sectors are given by Eqs. (jlV.lSj) . 

Case 5. A(z) = ^z(z 2 + 2(1 - 2m)z + l) . 

Change of variable: z = — — -. 

1 — cn q 

Supercharge: 



W 



cn q + 2 dn 2 q dn q 

h, = , t = -2 , 

sng dng sng 

(26i - N - 2a~ + 1) sn 2 q + 2{N + 2cC - 1] 



(IV.24) 



4 sng dng 
Potentials: 



T/± cr^a* ~ 1) ~ 1) A/" m/3 ± cn g + (/^) 2 + \ mm' {I - N 2 / 

V ~ 4(l + cng) + 4(1 -cng) + 2~d^q + V ° ' (IV - 25) 



Parameters: 

P ± = ±i(26i - (1 - 2m)(7V + 2a" - 1)) . (IV.26) 

The scalar potentials V ± are related by a rea/ translation, namely 

7+( g ) = V~(q + 2K). (IV.27) 

Hence if + and H~ are self-iso spectral in this case, cf. (2 
Gauge factors: 

u ± /l-cngx^ 1 / /3± Pcng + £; ,2 \ . , 

where k' = y/m/ = y/1 — m. 

The potentials V ± are singular at integer multiples of 2K, so that their corresponding 
Hamiltonians are naturally defined on, e.g., the interval S = (0, 2K). The necessary and 
sufficient conditions ensuring the square integrability of the wave functions in the solvable 
sectors are given by 

y± C L 2 (S) _I< a ±<Ar 2 -Ar 1 + ^, (IV.29a) 

V^ 2 CL 2 (S) <=^ - l -+M 2 -N 1 <a ± <^-. (IV.29b) 

Note, in particular, that the condition M\ > A2 implies that the inequalities (|IV.29ajl cannot 
hold unless N\ - N 2 = 0, 1. 
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V. AA-FOLD SUPERALGEBRA AND ASSOCIATED POLYNOMIAL FAMILIES 



In ordinary ("one-fold") supersymmetric quantum mechanics, it is well known that the 
superHamiltonian H is proportional to the anticommutator of the supercharges Q%-, and 
hence the operators and H span a three-dimensional Lie superalgebra. This cannot 
possibly be the case in TV-fold supersymmetric quantum mechanics (with JV > 1), since the 
anticommutator 

H*=±{Qfr,Qjr} (V.l) 

is now a linear differential operator of order 27V. Note, however, that from the nilpotency 
of the supercharges (Eq. (|II.5[) ) and the supersymmetric character of H (Eq. (jll.8|) h it 
immediately follows that the anticommutator (|V.1[) . that we shall henceforth call the mother 
Hamiltonian, commutes with Q^- and with H, namely 

[H N ,Q%\ = [Wa/-,H] =0. (V.2) 

These relations strongly suggest that the mother Hamiltonian TIm is a polynomial of degree 
/V in H, say TCj\f = IL\f(H), and hence that the operators Q^- and H span a nonlinear 
superalgebra of degree Af defined by the relations (jll.5|) . (|II.8jl . and 

{q^qj/} = 2u^(n). (v.3) 

That this is indeed the case was rigorously proved in Refs. P, HI , where it was also shown 
that ILv is proportional to the characteristic polynomial of the restriction of the component 

■V; 



Hamiltonians H ± to the invariant spaces Vtr, namely 



U M (E) = 2 N ~ l det(if ± | v ± - E) . (V.4) 
By Eqs. ()H.37j) and ()HI.35j) . this is equivalent to the "gauged" relation 

Hm{E) = 2 Ar - 1 det (H% ± - E) . (V.5) 



It was later shown in Ref. [1J] that for type A models the right-hand side of the latter equa- 
tion is proportional to the critical generalized Bender-Dunne polynomial (GBDP) associated 

to the gauged Hamiltonian H , which in this case is Lie- algebraic with respect to the Lie 
algebra sl(2). Bender-Dunne polynomials were introduced in Ref. to determine the 
solvablepart of the spectrum of a well-known quasi-exactly solvable sextic oscillator Hamil- 
tonian p, and were soon generalized by Finkel et al. [3i| to virtually all one-dimensional 
quasi-exactly solvable models associated to the sl(2) algebra. 

We shall now determine the polynomial for the type C models constructed in the 
previous section; in particular, we shall show that in this case Uj\f factorizes (up to a multi- 
plicative constant) as the product of two critical Bender-Dunne type polynomials of degrees 
M and Af 2 . In view of the equality of the characteristic polynomials of the operators H |y- 

andH + \y+ (which follows in general from Eq. (jV.5|) ). it suffices to study the action of one 
of the gauged Hamiltonians in its corresponding invariant space, e.g., of H = H~ in 

Recall, to begin with, that the linear space is the direct sum 

% = <1, z, ■ ■ ■ , ®z"(l,z,... } z^- 1 ) EE Vj£ } © < (V.6) 
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of the two spaces vffl (cf. Eq. (jIII.3)0 . each of which is invariant under the action of H. 
Hence 

det (H\v M - E) = det - E) ■ det (H^ - E) , (V.7) 

where we have set 

# W =% (A) . (V.8) 

The second-order operator ^~( l_1 ) A ^"W^(* -1 )^ (i = 1 ; 2) preserves a type A space (jIII.3|) of 
dimension Afi, and the coefficient of d 2 z in this operator is the polynomial A(z). Since A(z) 
has degree three and vanishes at the origin (see Eq. (|III.9ajl ). it follows that each of the 
operators z -^ x H^z^ x (i = 1,2) defines an associated family of GBDPs {vrj^ 1 '^ 2 '}^ 
whose critical element irfj}*^ 2 ^ (since Afi = dimVj^) is proportional to the characteristic 
polynomial of z~( l ~^ x H W^C*- 1 )^ and hence of H^: 

det(H (l) -E) = det(z- {i - 1)x H^z (l - 1)x - E) = (-l^vrj^ 2 ^) , % = 1, 2 . (V.9) 
By Eqs. (|V.5J) and ()V.7j) we thus have 

Iw = (-lJ 2 vr lM 7r 2M , (V.10) 
and consequently (cf. Eq. (|V.3jl ) the type C A/"-fold superalgebra is given by 

{Qir,Q%}= [<&,H] =0, (Vila) 
{Qat, Qjf} = (-2)^S^(H)vr^(H). (V.llb) 

We should note at this point that, as previously remarked, the operators of the form 
z - (i-i) \ z (i-i)\ ^ _ ^ 2) do not exhaust all the possible gauged Hamiltonians of type A. 
As a consequence, some of the characteristic features of the type A models are inevitably 
lost in the type C case. In particular, the GL(2,R) invariance of type A models, which 
ensures that all the coefficients of the associated GBDPs are expressed in terms of polyno- 
mial invariants !14], is broken in type C models and, as a consequence, the same is true for 
the Af-iold superalgebra of type C. In other words, the type C models have two invariant 
subspaces of type A at the cost of the GL(2,R) symmetry. On the other hand, it turns out 
that each of the polynomial families associated with type C models acquires a novel feature, 
namely the dependency on two positive integers Afi and A/2 (cf. Eq. (jV.18|) below). That is 
the reason for the rather cumbersome notation vr]"^ 1 "^ 2 ', that we shall hereafter abbreviate 
as 7Tj fc unless the dependence on Aft and A/2 is crucial. 

To construct the polynomial families {n'i l .}'kLo (i = 1,2) associated with type C models, 
let Xe{z) denote an eigenfunction of H with eigenvalue E. In view of Eq. (fV6jl . we shall 
consider the following two formal expansions of this eigenfunction in powers of z: 



00 



X e(z) = z'- l > x g fikTt+jhljX) Z " ' i=1 ' 2 ' <V ' 12) 



Here we have set itn- = Jik^ik, Jik being a numerical coefficient that must be chosen so 

)(A) 



that 7Tj fc is monic. Clearly, the necessary and sufficient condition for Xe( z ) to belong to vi A ^ 



is that 

7t ik (E) = , for all k > Afi . (V.13) 
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Acting on Xe{z) with H and using Eqs. (jII.14|) . ()III.9|) . and one immediately arrives 

at the following recursion relation for the coefficients TXi^{E): 



ax(k + 2) A + 1)^+1 = 

- E + c + (k + (i-l)X)(b 1 + a 2 (k + (i-l)X-M+l)) 



-a 3 (k + (i- 1)A) (k+(i- 1)A - l) (k + (i - 2)A - N + l) 



+ (M 1 -1)(M 2 + X-1) 



k > 



0. (V.14) 



Note that the coefficient a\ is nonzero in all the canonical forms listed in Table |l] with the 
exception of the second one, which corresponds to the trivial case of a constant potential. 
Hence we shall assume in what follows that a\ ^ 0. In that case the recurrence relation 
f|V.14j) can be brought to the more standard form 



E + c +(k + (i- 1)A) (6i + a 2 (k + (i - 1)A - AT + l)) 



ai a 3 (k+(i- 1)A) (A; + (i - 2) A) (ife + (z - 1)A - l) 



x (fe + (z-2)A-JV+l) + (M-l)(Wa + A-l) 



7Ti,fc— 1 j 



fc > 



(V.15) 



by choosing the so far undetermined multipliers 7^ as follows: 



7»,Jfc 



- ai ) fc r(A; + l + (2-2)A) 



fc > 0. 



(V.16) 



The three-term recursion relation ()V.15|) will actually define a family of weakly orthogonal 
polynomials {^i,k}^a provided that the coefficient of it^k-i vanishes for some non-negative 
integer value of k j3l|. If k = K is the lowest such value, and Ej (j = 1, . . . , K) is a root of 
the critical polynomia l 7ii k, then the recursion relation (jV.15jl implies that 7Ti t k(Ej) = for 
all k > K. It follows 29, 3CJ] that the linear space 1 - )A (1, z, . . . , z K l ) is invariant under 
H, and that the eigenvalues of the restriction of H to this space are the K roots (counting 
multiplicities) of the critical polynomial tt^k ■ in particular, the characteristic polynomial of 



the restriction of H to the invariant space z^ 1 - )A (1, z, 



K-l 



For the recursion relation (|V.15j) . the coefficient of Tti.k-i can be written as 

-ataskik - Mi){k - X)(k-N 2 - A), i = l 
- ai a 3 k(k - M 2 )(k + X)(k - Mi + X) , i = 2. 



) is proportional to tz^k{E). 



(V.17) 



Taking into account the restrictions (|III.28|) on A, it is easily seen that the degree of 
the critical polynomial is K = Mi, i = 1,2. This establishes Eq. (IV. 9)1 . since the space 



z (% 1)^(1 ; z , . . . , z Mi 1 ) coincides with V 



(A) 

Mi 



1,2) by Eq. (jVTfijl . 



It is worth mentioning that the polynomial systems n^f. associated with type C A/"-fold 
supersymmetry are always weakly orthogonal in spite of the fact that the solvable sectors 
Vjy- are not always normalizable, as was studied in the preceding section. We thus obtain 
further evidence of the claim in Ref. jlj] that normalizability has in general nothing to do 
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with the weak orthogonality of the associated GBDPs. Note also that Eq. (jV.15|) becomes a 
two-term recursion relation if and only if the coefficient 03 vanishes, or, equivalently, if the 
corresponding Hamiltonian is solvable in Turbiner's sense. 

The recursion relation (|V.15J) determining the Bender-Dunne type polynomials n^k and, 
ultimately, the TV-fold superalgebra via Eq. (jV.llbj) . can be recast into the following more 
concise form using Eq. ()V.17j) : 



E + c + (k + (i - 1)A) (61 + a 2 (k + (i - 1)A -Af+ l)) 



ai a 3 k(k - Mi) (k + (2i - 3) A) (k - A/" 3 -i + (2i - 3)A)7r ijfc _ 1 . 



(V.18) 



This recursion relation can be used without difficulty to compute the polynomial families 
^Ffc 1 anc ^' ^ n particular, the critical polynomials determining the A/"-fold superalgebra, for 
any given values of Al and M 2 - We shall exhibit in what follows a few examples of these 
polynomials for M = 3 and 4, assuming that A/1 and A/2 satisfy the restriction A/1 > A/2 
imposed in Section ITTT1 

Example 1. A/1 = 2, A/" 2 = 1. 
Polynomial system: 

nff(E) 

7i 



(V.19a) 
(V.19b) 

(V.19c) 



E + c Q , 

(E + cq)(E + c + 61 - a 2 ) + A(A - l)aia 3 
E + c + A6i + A(A-2)a 2 . 

3- fold superalgebra: 

{Qt,Qs} = -8[(H + co)(H + c + 6i-a2) + A(A-l)a 1 a 3 ] [H + c + A&i + A(A-2)a 2 ] . (V.20) 

Example 2. A/1 = 3, N 2 = 1. 

Polynomial system: 

7iff(E) = E + c , (V.21a) 

4 3 2 1] (^) = + co)(-E + c + 61 - 2a 2 ) + 2A(A - l)aaa 3 , (V.21b) 
7r5 1] (^) = (# + c )(^ + c + 61 - 2a 2 )(E + c + 2b x - 2a 2 ) 

+ 4(A - l)aia 3 ((A - 1)(E + c ) + A(6i - a 2 )) (V.21c) 

7r^i 1] (S) =E + c + A61 + A(A - 3)a 2 . (V.21d) 

4- fold superalgebra: 

{Qt> Qa } 
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H + c )(H + c + 61 - 2a 2 )(H + c + 2b x - 2a 2 ) 

+ 4(A - l) Ql a 3 ((A - 1)(H + c ) + X(h - a 2 )) 
x [H + c + A&i + A(A - 3)a 2 ] . 

Example 3. A/1 = 2, A/ a = 2. 



(V.22) 
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Polynomial system: 



n [ ?f(E)=E + c 



IT 



[2,2] 
1,2 



(E) = (E + c ){E + c + 61 - 2a 2 ) + (A 2 - l)a x a 3 . 



4 2 f ] (^) = # + co + A61 + A(A - 3)a 2 . 
Tr l ?f(E) = (E + c + (X + l)h + (A + 1)(A - 2)o 2 ) 

x (E + c + A61 + A(A - 3)a 2 ) + (A 2 - l)aio 3 



4-fold superalgebra: 

{Qt,Ql} 



16 [(H + c )(H + c + 61 - 2a 2 ) + (A 2 - 1)0^3] 
x (H + c + (A + l)6 1 + (A + l)(A-2)a 2 ) 

x (H + c + A61 + A(A - 3)a 2 ) + (A 2 - 1)0103 



(V.23a) 
(V.23b) 

(V.23c) 
(V.23d) 



(V.24) 



VI. SUMMARY AND DISCUSSION 

In this article we develop an algorithmic procedure for constructing an A/"-fold supersym- 
metric quantum system starting from a given finite-dimensional space of functions invariant 
under a suitable gauge transform of one of the supercharges. Although the method is very 
general, it is especially useful when one knows a particular example of quasi-solvable Hamil- 
tonian and its algebraically solvable wave functions. We have applied this procedure to the 
monomial spaces of Post-Turbiner type, thus obtaining the new models of type C as well as 
recovering the previously known type A and type B systems. We would also like to stress 
that, although the procedure always yields an A/"-fold supersymmetric system {H ± 1 P J ^) 1 it 
does not rule out the existence of more general systems of the same type. For instance, 
for the type C models specifically discussed in this paper the supercharges are given by 
Eqs. (III.3J) and (|III.22J) . and the method developed in this paper guarantees that the su- 
persymmetry algebra holds for suitable Hamiltonians provided that the functions E 
and F satisfy Eq. (|III.18J) . It is not clear, however, whether this sufficient condition is also 
necessary. It should also be noted in this respect that the situation is completely analogous 
for the type B models discussed in Ref. |l5| . 

The normalizability of the solvable sectors of the type C models, which plays an impor- 
tant role for the existence of dynamical A/"- fold supersymmetry breaking [l|, [l2|, is briefly 
investigated in Section |IV| From Eqs. (jTV^jl . (llV.15j) . (llV.18j) . and 1I1V.29|) . we see that 
there is little chance for both sectors Vjy- and Vjjr (resp. Vjj- and Vjy ) to be simultaneously 
normalizable. This means that, in general, only a part of the whole solvable sector Vtf 
(resp. Vjj) is physical in type C A/"-fold supersymmetric systems. For example, among the 
four subsectors V^. in Case 1 only Vjj is normalizable and thus physical (boundary condi- 
tions aside) if b\ < and A > 1. We can thus say that A/"-fold supersymmetry is partially 
broken in type C models. This phenomenon is novel for type C due to the characteristic 
structure of its solvable sector, and, to the best of the authors' knowledge, has not been pre- 
viously mentioned in the literature. For a more precise discussion, we need of course to deal 
with the boundary conditions at the singularities in order to define a self-adjoint extension 
of the Hamiltonians. This kind of mathematical subtlety is beyond the scope of the present 
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article, and we shall therefore content ourselves with referring the reader to, e.g, Ref. [32j 
for a recent discussion of this topic. We also note, in this connection, that the significance of 
the boundary conditions in (ordinary) supersymmetry breaking was recently reported in a 
different context, namely, through the careful calculation of the fermion determinant arising 
out of the path integral formalism 33] . 

The structure of the solvable sectors in type C models gives rise to another interesting 
phenomenon in the theory of exact solutions of the Schrodinger equation. Indeed, when 
the parameter a + (resp. or) is a positive integer I + 1 (in both cases A is a half-integer, 
cf. Eq. (IIII.33|) ). the potential V + (resp. V~) in Case 1 given by Eq. (jIV.4|) is nothing but the 
radial harmonic oscillator plus the centrifugal potential with a properly quantized angular 
momentum I. It is then apparent that for b\ > (resp. h\ < 0) the first solvable sector 
Vj} (resp. Vjj- ) (Mi = 1, 2, . . . ) in Eq. (|IV.5ajl contains the physical solutions around q = 0, 
while the second sector Vjt (resp. Vjj- ) in Eq. (jIV.5b|) yields the "second" solutions, usually 
discarded as unphysical due to their singularity 5 at the origin. 

By Eq. ()HI.36j) . we may in fact say that type C models are characterized by the fact 
that both linearly independent solutions around z(q) = are quasi-solvable. More precisely, 
from Eqs. (pT45j) . ill I I.Qcj) . ill 1.1.1 0)1 . (jl ll.lfij) . and Table [fl it is easily seen that the operator 
Hjj- — E (resp. Hj^- — E) has a regular singularity at z — 0, so that Fuchs's theorem applies. 
Except in the trivial Case 2, which we shall henceforth ignore, for all values of the energy 
E the roots of the indicial equation are and A (resp. Mi and Mi + A). From Eqs. (jIII.13j) 
for the space Vjf, and the analogous decomposition 

= <1, . . . , z*~ l ) © z^ + ~ x (1, . . . , z*- 1 ) 

for V^, we see that each of the four sectors of the invariant spaces Vj^ consists of the solvable 

eigenfunctions of behaving as z l near z = 0, where i is one of the roots of the indicial 
equation listed above. Note, in this respect, that this remark is still valid when the difference 
of the roots of the indicial equation, which is given by A (resp. A), is an integer, even if in this 
case Fuchs's theorem can only guarantee the existence of one linearly independent power 
series solution around z = 0. 

The construction of the Bender-Dunne type polynomial systems associated with the 
type C models turned out to be straightforward. The breakdown of the GL(2, M) symmetry 
in type C models spoils the characteristic feature possessed by the type A polynomials. 
Instead, the GBDPs of type C have a novel dependence on two positive-integer parameters 
Mi and Mi- For any given pair (Mi, Mz) G N x N, we obtain two related families of weakly 
orthogonal polynomials satisfying the recursion relation ()V.18[) . The assertion in Ref. [1J| 
that normalizability has nothing to do with the weak orthogonality of the associated family 
of polynomials has also been confirmed. We also note that a similar construction for the 
type B systems presents considerable difficulties due to the lack of form-invariance under 
projective transformations of these models. This difficulty has also prevented the systematic 
calculation of the explicit form of the type B A/"-fold superalgebra. 

All the quasi- solvable second-order differential operators in one variable preserving a 
finite-dimensional linear space of monomials, which were classified in Ref. [lfij . have now 
been brought into the framework of A/"-fold supersymmetry. 6 The natural continuation of 



5 Or their non- vanishing, for I = 0. 

6 With only a few exceptions for TV = 3 and TV = 4. 
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the present work is the study of operators possessing non-monomial type invariant sub- 
spaces, and the construction of the associated supersymmetric models following the general 
algorithm developed in Section |TlJ It would also be of great interest to extend the results 
obtained in this paper to a multi-dimensional space-time. Indeed, the generalization of M- 
fold supersymmetry to higher-dimensional space-times remains one of the most challenging 
open problems in this field. Recently, Smilga showed jjjij that some A/"-fold supersymmetric 
systems with M = 2 can be realized as weakly supersymmetric field theories in one space- 
time dimension. Although the results in Ref. |34j clearly indicate the existence of rather 
severe obstructions, further investigation in this direction would still be certainly worth 
undertaking. 
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